Observation of an ABC effect in proton-proton collisions 
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The cross section for inclusive multipion production in the pp — > ppX reaction was measured 
at COSY-ANKE at four beam energies, 0.8, 1.1, 1.4, and 2.0 GeV, for low excitation energy in 
the final pp system, such that the diproton quasi-particle is in the So state. At the three higher 
energies the missing mass Mx spectra show a strong enhancement at low Mx , corresponding to an 
ABC effect that moves steadily to larger values as the energy is increased. Despite the missing-mass 
structure looking very different at 0.8 GeV, the variation with Mx and beam energy are consistent 
with two-pion production being mediated through the excitation of two A(1232) isobars, coupled 
to S- and D— states of the initial pp system. 

PACS numbers: 13.60.Le 13.75.Cs 14.40.Aq 25.40.Qa 



The ABC effect is a sharp enhancement of the two- 
pion invariant mass spectrum near its threshold that has 
puzzled hadron physicists for almost fifty years. First 
observed in a pd — > 3 HeY missing- mass experiment 1|, 
and subsequently confirmed using a deuteron beam [2|, 
it was seen most spectacularly in the dd — > 4 HeA reac- 
tion 0,13. The ABC typically manifests itself as a peak 
at an invariant mass Mx ~ 310 — 330 MeV/c 2 and it was 
first even speculated that it might be a scalar meson with 
isospin I — 0. However, the fact that the mass and width 
varied with experimental conditions [2j suggests that the 
ABC is a kinematic enhancement. This is consistent with 
the smooth behavior of the isoscalar tttt phase shifts at 
low energies. To understand the nature of the ABC it is 
necessary to investigate the effect in simpler systems. 

The ABC showed up very clearly in the zero-degree 
momentum spectrum of the deuteron from np — > dX [B[ . 
Here it leads to two peaks, i.e. forward and backward pro- 
duction in the cm. system, though the momentum spread 
of the neutron beam degraded the missing-mass resolu- 
tion. This has recently been overcome by measuring the 
quasi- free pd — » d7r 7r°p sp reaction semi-exclusively at 
beam energies of T p = 1.03 and 1.35 GeV, with just the 
spectator proton p sp escaping detection Q. The only 
published study of the ABC effect in proton-proton col- 
lisions was carried out at T p = 1.52 and 1.81 GeV Q, 
though this concentrated on possible substructure. 

Studies of the energy dependence of the pd — > 3 HeY 
and dd — > 4 HeY cross sections showed that the ABC 



effect is most prominent for energies where the maxi- 
mum missing mass is about 600 MeV/c 2 @|. Since this 
corresponds to the mass difference between two A(1232) 
isobars and two nucleons, it is tempting to suppose that 
the two-pion production is mediated by the excitation 
and decay of two separate A isobars. Such a model with 
7T Q and then n + p exchange Q had semi-quantitative 
success in describing the existing np — > dX data. The 
dominant contribution to the cross section arises when 
the two A are in a relative 5-wave. For the production 
of an /tttt = pion pair, isospin conservation requires 
Iaa = 0. It follows from the generalized Pauli princi- 
ple that the two isobars must then couple to a total spin 
of Saa = 1 or 3. In contrast, for pp — > AA — > ppirir, 
Iaa = 1 so that Saa = or 2, and an s-wave pion pair 
can have 7 T7r = or 2. The investigation of the ABC 
phenomenon in np and pp collisions is therefore largely 
complementary. 

The aim of the present experiment is the measurement 
of the pp — > {pp} s X reaction when the excitation energy 
in the final diproton system is low, E pp < 3 MeV. This 
quasi-particle, denoted by {pp} s , will be almost exclu- 
sively in the ^So state. Under these conditions, the kine- 
matics are similar to those of pn — > dX, which facilitates 
a comparison of the two reactions. Of the four proton 
beam energies used, T p — 0.8, 1.1, 1.4, and 2.0 GeV, the 
lowest is well below the nominal AA threshold, and by 
2 GeV much of the A strength has passed. 

The measurement of single and multipion production 
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in pp collisions was performed using the ANKE mag- 
netic spectrometer [HJ, positioned at an internal tar- 
get station of the COSY synchrotron storage ring of the 
Forschungszentrum Jiilich. The single pion results have 
already been reported The proton beam was di- 

rected at a hydrogen cluster-jet target with an areal den- 
sity of 2 x 10 14 atoms/cm 2 . The resulting positively 
charged particles were recorded in the ANKE forward 
detector system, consisting of three MWPCs and a two- 
layer scintillation hodoscope, the whole giving a horizon- 
tal acceptance of 12° and ±3.5° vertically. The parti- 
cle momenta, deduced by tracing through the analyzing 
magnetic field, had typical uncertainties of w 1%. The 
emerging proton pairs were identified by evaluating the 
difference At between the times of flight measured with 
the hodoscope and those calculated from the particle mo- 
menta, assuming the particles to be protons [111 ]. 



tion shows a rapid rise from the two-pion threshold with 
a gentle decrease towards the upper kinematic limit. 




FIG. 1: Distribution in the square of the missing mass of the 
pp —* ppX reaction at 1.4 GeV. Imposing a At cut to se- 
lect the two protons gives the open histogram. The lightly 
shaded plot corresponds to the additional requirement that 
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< 3 MeV and the heavy shading reflects the further 
'pp > 0.95 cut to match the overall ANKE acceptance. 
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The positions of the n° and r\ peak are indicated, as is the 
two-pion threshold (dotted line). 

Figure [T] shows the spectrum in the square of the miss- 
ing mass of the pp — > ppX reaction at 1.4 GeV. Selecting 
pp events on the basis of the At cut eliminates almost all 
d,Tr + pairs, leaving a background of less than 0.1%. To en- 
sure that the two protons are dominantly in the So state, 
a further cut on the pp excitation energy, E pp < 3 MeV, 
was imposed. After weighting by detection efficiency, 
such pairs are distributed isotropically in their rest frame, 
as expected for an 5-wave. The E pp distributions also The ABC position moves steadily to the right and gets 



FIG. 2: Distribution in missing-mass squared for the pp — ► 
{pp} s X reaction for E pp < 3 MeV and cos $ pp > 0.95 at a) 
0.8, b) 1.1, c) 1.4, and d) 2.0 GeV. The <q signal is seen at 
the expected position for the two higher energies. The curves 
represent normalized simulations within a phase-space model. 



The development of the counting rate distribution with 
energy is displayed in Fig. [5] for the multipion region. 
Also shown there are the results of a phase-space Monte- 
Carlo simulation of the detector and analysis code, in- 
cluding effects of momentum resolution, and applying 
the same cuts as for the experiment. In no case does 
this approach reproduce the data. 

A simulation was used to deduce the differential cross 
sections of Fig. [3] from the counting rate data, taking as 
input the double- A model to be discussed with Eq. ([3]). 
The simulation also defines the momentum reconstruc- 
tion uncertainties and this leads to a ir° peak whose width 
is compatible with that shown in Fig. [T] The requisite 
luminosity was derived from the simultaneous measure- 
ment of pp elastic scattering 11 1 . 

As shown in Fig. [31 the strength at 0.8 GeV is pushed 
towards the maximum missing mass, while at the higher 
three energies there is a peak of the ABC type at low 
Mx and no sign of any enhancement at large masses. 
Above its threshold there are also contributions from rj 
production, with cross sections integrated over our cuts 
of a„ = 4.3 ± 0.8 at 1.4 GeV and 4.5 ± 0.9 nb at 2.0 GeV. 



follow closely the expected 5-wave final-state-interaction 
behavior. Since for low missing masses the cm. angle 
of the diproton, $ pp , is strongly limited by the ANKE 
acceptance, we imposed the extra cut cosi? pp > 0.95 on 
the whole spectrum. Apart from peaks corresponding to 
the production of single 7r° and rj mesons, the distribu- 



broader as T p increases above 1.1 GeV, a confirmation 
that the phenomenon is primarily a kinematic effect. The 
peak seems to be rather clearer than that observed at 
Saclay at 1.5 and 1.8 GeV Q. Although a previous ex- 
clusive measurement of pp — ► ppir + T:~ at 0.8 GeV [l3| 
did not study in detail our kinematic configuration, the 
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FIG. 3: The pp — > ppX differential cross section with statis- 
tical errors as a function of the square of the missing mass at 
a) 0.8, b) 1.1, c) 1.4, and d) 2.0 GeV for E pp < 3 MeV and 
cosi9 pp > 0.95. The 77 peaks are indicated. Events were sim- 
ulated using Eq. @ with Ad = (long dashes), with As = 
(short-dashed), and with the best fit of Table U (solid line). 



two data sets are consistent in the region of overlap. 

The total cross section for the production of three pi- 
ons in pp collisions at 1.35 GeV is about 0.6% of that 
for two pions [13] and so it is only at 2 GeV that there 
might be any significant contribution to the data from 
the pp — > ppmnr reaction. The modelling of the results 
must therefore be based upon the assumption that two- 
pion production dominates. 

A AA model applied to two-meson production in the 
pn — > dX or dd —> aX reaction allows the momentum 
transfer to be shared between the nucleons, facilitating 
the formation of the observed deuteron or a-particle. The 
strong peaking of the A decay in the forward and back- 
ward directions also enhances this sticking factor. It is 
therefore particularly suited to our conditions but might 
not be appropriate for other kinematic domains. A full 
estimation of two-pion production in such a model de- 
pends sensitively upon the exchange mesons and form 
factors etc. Nevertheless, the general structure can be 
written down if we assume that the intermediate A are 
in a relative S-wave so that the initial protons are in the 
^So or 1 Z?2 wave. Since the p-wave nature of each A de- 
cay necessarily involves a momentum factor, the obvious 
ansatz for the pp — > {pp} s irir matrix element is 



M 



A s (k[ ■ k' 2 ) + \A D (3(p • k[)(p ■ k' 2 ) - k[ ■ k' 2 ) , (1) 

where p is the cm. proton beam momentum and k[ is the 
relative momentum of the i'th pion with respect to one of 



the recoiling protons. The S and -D-wave amplitudes As 
and Ad are scalar functions that may depend strongly 
upon the kinematic variables because of the excitation of 
the two A, though we shall neglect any variation with 
Mx- In terms of the relative irn momentum, 2q, and 
that of the diproton, —k, the pion cm. momenta are 
ki_2 = \k ± q, and the matrix element becomes 



M = A s (a 2 k 2 - (3 2 q 2 ) 
+ \A D (3(a 2 k 2 z -(3 2 q 2 z )~(a 2 k 2 -(3 2 q 2 )), 



(2) 



where the z-direction is taken along p. In this approxi- 
mation, q 2 and k 2 are linked by energy conservation. For 
simplicity of presentation, a Galilean transformation is 
used to evaluate the recoil factors a = \{m + 2^) / {m + p) 
and (5 = m/(m + /1) in terms of the pion (/i) and proton 
(m) masses. Only the part of the matrix element pro- 
portional to 3q 2 — q 2 corresponds to a c£-wave in the 7T7r 
system and this contributes very little to the ABC peak, 
which is s-wave in nature. 

Since the direction of the 7T7t relative momentum is 
not observed in a missing- mass experiment, \A4\ 2 must 
be averaged over the angles of q; 



2;„2 



< \M\ 2 > = \A s (a 2 k 



j3 2 q 2 ) + \A D a 2 {?,k 2 z ~k 2 )\ 



(3) 



Figure [3] shows the results of implementing the simple 
double- A model of Eq. ^ . Since the charges of the pi- 
ons were not measured, we assumed an 7 W7r = ratio of 
7r + 7r~ : 7r°7r° = 2:1, with the corresponding masses being 
used in the two phase spaces. If there were only S'-wave 
production, < \A4\ 2 > would be maximal at q 2 — (the 
ABC bump) and at k 2 = (the maximum allowed miss- 
ing mass), with a deep valley between these two struc- 
tures. Such behavior is seen for dd -> 4 HeX @,||. The 
Ad = case is illustrated in Fig.[3j as is that of a pure D- 
wave, which gives a broad maximum at low 7T7t invariant 
masses and no sign of any high mass bump. 

The general shapes of the spectra at the three higher 
energies are qualitatively reproduced by the pure D-wave 
model. However, the peaking towards maximum Mx in 
the 0.8 GeV data shows the necessity for a substantial 
S'-wave contribution at this energy. To investigate this, 
the uncorrected data of Fig. [5] were fitted to the Monte- 
Carlo simulations using as input Eq. §S§ with free values 
of As and Ad- The values of these parameters that best 
reproduce the spectra are given in Table U along with the 
resulting integrated cross sections. Only at 0.8 GeV is 
As/ Ad required to be large. The description of the data 
shown in Fig. [3]is much improved, especially at 0.8 GeV, 
but the model does not reproduce the sharpness of the 
ABC peaks seen at 1.1 and 1.4 GeV. The values of Ad fall 
steadily with energy but, to understand the significance 
of this, would require a full dynamical model. 

The special quasi-particle kinematics studied here are 
very similar to those of pn — > d7r°7r° and so a com- 
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TABLE I: Cross section ax for the pp — > {pp} s X reaction 
integrated over cos$ pp > 0.95, E vv < 3 MeV and over Mx 
up to the kinematical limit. The experimental acceptance is 
limited to less than M™ ax , and ax is obtained by correcting 
the measured cross-section by an amount £, as suggested by 
fits of the AA model to the differential data. In addition to 
the given statistical errors, there are systematic uncertainties 
of 7%, coming principally from the luminosity and acceptance 
evaluations [f]J |. The parameters As and Ad are determined 
by fitting the spectra away from regions of rapid variation, 
where resolution questions arise, and at too high masses to 
avoid the tail from p production. The results for Ad are 
normalized at 1.1 GeV. 



T v 


M m ax 


ax 




A s /A D 


A D 


GeV 


MeV/c 2 


nb 


% 






0.8 


351 


12 ±2 


34 


-1.23 ± 0.10 


1.09 ±0.05 


1.1 


468 


389 ±6 


2.3 


-0.20 ±0.03 


1.00 


1.4 


571 


563 ±6 


2.4 


-0.10 ±0.03 


0.48 ±0.01 


2.0 


741 


456 ±5 


5.8 


-0.23 ± 0.02 


0.17 ±0.01 



parison with the CELSIUS- WASA data 6] is instruc- 
tive because of the different spin-isospin selection rules. 
The measurement of the Fermi momentum when using 
a deuterium target allowed the determination of cross 
sections as a function of the laboratory kinetic energy 
T p . The results suggested a resonance-like structure with 
a mass of 2.39 GeV/c 2 and width of T « 90 MeV/c 2 
(T p w 1.17 GeV, T T = 0.23 GeV). Though our energy 
intervals are rather wide, no such behavior is obvious in 
the values given in Table Q] which would suggest a strong 
isospin dependence of the production mechanism. 

Within the S-w&ve AA approach, the matrix element 
of Eq. {l} allows for both As and Ad terms. If two- 
pion production were driven instead by the excitation of 
5-wave Y*(1440)Y pairs, where N* is the Roper reso- 
nance, an As term is generated through the double p- 
wave cascade decay Y*(1440) — > A7r — > pirir. However, 
this model requires that Ad — 0. Contributions from 
both N*N and A A might therefore provide an expla- 
nation for the energy dependence of As /Ad shown in 
Table fl] It should be noted that the AA model would 
lead to an ABC effect also in pp — > {nn} s Tr + ir + , which is 
forbidden in the N*N model. Earlier studies of the ABC 
effect were not sensitive to any I n „ — 2 component. 

In summary, we have measured the differential cross 
section for the pp — > {pp} s X reaction at four beam en- 
ergies from 0.8 to 2.0 GeV under the specific kinematic 
conditions where the proton-proton excitation energy is 
below 3 MeV and the cm. angle between the diproton 
momentum and the beam axis is less than 18°. The form 
of the Epp spectra and isotropy of the angular distribu- 
tion in the diproton frame are consistent with the two 
final protons being in the ^So state. Strong deviations 
from phase space are observed in the missing-mass vari- 



able, with a peak in Mx whose position varies with beam 
energy, and no broad bump at maximum missing mass of 
the type observed in dd — > aX 0, 0] • On the other hand, 
there is no sign of the resonance-like structure in the en- 
ergy dependence of the integrated cross section that is 
reported for pn — > d7r°7r° [6|. 

The evidence presented here points towards to the 
ABC produced in pp collisions being a kinematic effect 
and this is backed up by our naive phenomenological de- 
scription. The isospin structure of pp — > {pp} s X is much 
richer than in the pn — > dX case and it is quite likely that 
there will be significant contributions from both I nn = 
and 2; a kinematic enhancement does not have to have 
a definite isospin. On the other hand, the spin structure 
is much simpler and, after assuming that each of the two 
pions is emitted in a p-wave, there are only two ampli- 
tudes to be considered. Except for the 0.8 GeV data, 
the bulk of the Mx spectra are consistent with a domi- 
nant D-wave coupling to the incident protons, whereas 
at the lowest energy the S'-wave is significant. The ma- 
jor discrepancies in the model are in the heights of the 
ABC peaks at 1.1 and 1.4 GeV, but a fully microscopic 
approach would be required to assess the origin of these. 

Exclusive measurements of pp — > {pp} s (irir) a over a 
wider range of angles would provide more stringent tests 
of the simple phenomenological description and the study 
of pp — > {rm} s 7r + 7r + would be particularly valuable. 
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